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We describe a simple magneto-optical experiment and introduce a magneto-optical Kramers- 
Kronig analysis (MOKKA) that together allow extracting the complex dielectric function for left- 
and right-handed circular polarizations in a broad range of frequencies without actually generating 
circularly polarized light. The experiment consists of measuring reflectivity and Kerr rotation, or 
alternatively transmission and Faraday rotation, at normal incidence using only standard broadband 
polarizers without retarders or quarter-wave plates. In a common case, where the magneto-optical 
rotation is small (below ~ 0.2 rad), a fast measurement protocol can be realized, where the polarizers 
are fixed at 45° with respect to each other. Apart from the time-effectiveness, the advantage of this 
protocol is that it can be implemented at ultra-high magnetic fields and in other situations, where an 
in-situ polarizer rotation is difficult. Overall, the proposed technique can be regarded as a magneto¬ 
optical generalization of the conventional Kramers-Kronig analysis of reflectivity on bulk samples 
and the Kramers-Kronig constrained variational analysis of more complex types of spectral data. 
We demonstrate the application of this method to the textbook semimetals bismuth and graphite 
and also use it to obtain handedness-resolved magneto-absorption spectra of graphene on SiC. 

PACS numbers: 78.20.-e,78.20.Ci,78.20.Ls,78.40.Kc,78.67.Wj 


I. INTRODUCTION 

Applying a magnetic field to materials influences pro¬ 
foundly their optical properties by the Zeeman splitting, 
formation of Landau levels (LLs), modification of the op¬ 
tical selection rules and sometimes inducing phase tran¬ 
sitions. The strength of magneto-optical (MO) effects 
depends on many factors, such as the density and mo¬ 
bility of the charge carriers, the value of magnetic mo¬ 
ments and their interaction, the band structure and the 
spin-orbit coupling. Magnetic field breaks the time rever¬ 
sal symmetry (unless it is already broken spontaneously) 
and triggers non-reciprocal optical phenomena, such as 
the Faraday 1 and the Ken® rotation. The MO rotation, 
which was historically the first compelling evidence of the 
electromagnetic nature of light, plays nowadays a crucial 
role in telecommunications, data storage, laser technol¬ 
ogy and material characterization. 

It is well-known that magneto-resistivity and Hall- 
effect measurements are in general more informative than 
the DC zero-field transport. Likewise, the magneto¬ 
optical response carries more ample information about 
a physical system than zero-field optics. It is also more 
complicated from the experimental and theoretical point 
of view. While the linear electromagnetic response of an 
isotropic material without field is described by a scalar 
function, such as the complex optical conductivity cr(ui) 
or the dielectric function e(w) = 1 + 4iria(oj)/u], its MO 
counterpart is essentially a 3 x 3 tensor. The Onsager 
reciprocity relation^ impose on it the following struc¬ 
ture: 


instead of one (the field B is along the z axis). If the 
light polarization is confined in the xy plane, which is 
the case for normal-incidence experiments in the Faraday 
geometry, where the held is parallel to the propagation of 
light), then only the two components, e xx and e xyi play a 
role. This allows us to restrict our discussion to the 2x2 
tensor (we omit the B argument for brevity): 


Cxyifo') 

^xx(^) 


( 2 ) 


which has the eigenvalues e±(w) = e xx (u >) ± ie xy {u>) cor¬ 
responding to the right- and the left-handed circular po¬ 
larizations. It is notable that this reduction is valid not 
only for fully isotropic but also for uniaxial media with 
the optical axis along the z direction, which increases 
dramatically the number of materials to which our dis¬ 
cussion applies. 

The goal of a complete linear-response optical- 
spectroscopy experiment in zero field is to determine the 
real and the imaginary parts of the dielectric function. 
While certain techniques, such as spectroscopic cllipsom- 
etry or time-domain spectroscopy, can provide experi¬ 
mental access to both of them, other methods rely on 
the Kramers-Kronig (KK) relation^ ■ due to a limited 
set of measured optical quantities. A widespread tech¬ 
nique of this kind is reflectance spectroscopy, where one 
measures the absolute reflectivity R(u>) in a broad range 
of frequencies and determines the phase 9(ui) of the com¬ 
plex reflectivity coefficienlP. 


e(w,H) 


e. X y(uj,B^ 0 

~ e xy (w,B) ^XX (w,B) 0 

0 0 c zz (uj, B) 


( 1 ) 


r{oj) = 


1 - y/eM 

1 + -\/ e M 


( 3 ) 


which involves three complex functions, e xx , e xy and e zz , 


Jahoda found^ that the KK relation can be applied to 
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the function f(u) = lnr(w) = In y 7 R(u) + iO(uj): 

n . . 1 In \J~R{x) 

6(uj) = - P / - -dx + 7T, (4) 

n J~oo a; - w 

which can be used to restore the unknown complex phase 
by a direct integration of the measured reflectivity prop¬ 
erly extrapolated outside the experimental range. Due to 
the parity relation R(—oj) = R(w) this formula reduces 
to 


6{u) 



lni?(a;) 
x 2 — w 2 


dx + 7r. 


( 5 ) 


The knowledge of the phase together with the reflectivity 
itself allows, via the Fresnel equations or other boundary- 
conditions, the determination of other optical functions, 
including the dielectric function and the optical conduc¬ 
tivity. 

In magneto-optics, a complete experiment should al¬ 
low determination of the real and the imaginary parts 
for both e xx (ui) and e xy (w) or equivalently for both 
e±(w) = 1+47T ia±(u>)/ui. Obviously, the number of quan¬ 
tities to be measured should be doubled as compared to 
the case of B = 0. At a first glance, one should simply 
work with circular polarized light and apply the zero-field 
techniques to each handedness separately. Although con¬ 
ceptually appealing, this strategy is difficult to realize in 
broadband spectroscopy as it is experimentally challeng¬ 
ing to generate circular polarized light in a broad range 
of frequencies. Standard waveplates, typically made of 
quartz, operate at selected wavelengths. There are re¬ 
ports of multilayer waveplates that cover a continuous 
range of f requencies, which is however still somewhat 
limitect-dih-4 Many waveplates must therefore be com¬ 
bined to cover a reasonably broad spectrum of circular 
polarized light, which is especially hard in the range of 
the phonon absorption in quartz. An alternative solution 
using elliptically polarized light was also demonstrated^ 
that allows reducing the number of waveplates but sig¬ 
nificantly complicates the analysis. A broadband phase- 
shifting element (retarder) using an internal reflection in 
silicon and other materials was fabricated 13 ^. Working 
with it is however associated with extra light absorption, 
the need of a delicate alignment and the penalty (in the 
case of Si) of not seeing visible light after the retarder. 

An alternative broadband solution without generating 
circular light is based on polarimetry, where one can mea¬ 
sure the MO rotation and ellipticity, or related quanti¬ 
ties. A rather sensitive polarimetric method is based on a 
photoelastic modulation (PEM), where one can measure 
these quantities using a lock-in demodulation at the main 
PEM frequency and at its second harmonic 15 15 1 How¬ 
ever, in order to achieve an experiment that is complete in 
the sense indicated above, one has to measure addition¬ 
ally the complex reflectivity or transmission. Further¬ 
more, the calibration procedure for the PEM-based mea¬ 
surement is quite involved. Another relevant technique is 
a spectroscopic magneto-ellipsometry with an access to 


many Muller-matrix elements 17 . While allowing a num¬ 
ber of independent quantities to be measured separately, 
this experiment unavoidably mixes the xy and z compo¬ 
nents of the 3x3 MO tensor due to the finite angle of 
incidence and thus adds extra complexity to the analysis. 
To our knowledge, a model-independent extraction of the 
complex functions e+(w) and e_(w) by this method has 
not been demonstrated yet. 

In this paper we propose a rather simple and fast 
pathway towards a complete broadband magneto-optical 
spectroscopy of materials with the MO conductivity de¬ 
scribed by Equation ([2]). It involves the measurement of 
the reflectivity (or transmission) and the Kerr (or Fara¬ 
day) rotation at a normal, or near-normal angle of in¬ 
cidence using two conventional polarizers without the 
need to measure other quantities, such as the elliptic¬ 
ity. The treatment of this set of data relies on a general¬ 
ized magneto-optical KK analysis (MOKKA), which we 
introduce and mathematically justify in this paper. 

To emphasize the basic ideas of the method, we shall 
restrict our consideration to materials, where the mag¬ 
netic and magneto-electric susceptibility are negligible as 
compared to the dielectric function, which is valid in a 
vast majority of non-magnetic compounds. In the dis¬ 
cussion, we shall however speculate on a possible inclu¬ 
sion of these electromagnetic terms that play an essential 
role in magnetic material^®, multiferroicd 1 ^ and topo¬ 
logical insulators^. For the sake of simplicity we shall 
also assume that the Fresnel equations are satisfied, even 
though the essential part of MOKKA holds even in the 
non-local limit, where more general electromagnetic re¬ 
lations should be used. 

The rest of the paper is organized as follows. Section 
|II| is devoted to the generalization of the conventional 
Kramers-Kronig reflectivity formalism to the case, where 
the magnetic field is set perpendicular to the sample 
surface. In section IIIIl the variational KK-constrained 


describes the experimental technique. Three application 
examples are presented in section [Vj The results are sum¬ 
marized in Section [VII 


analysis^ is extended to MO experiments. Section IV 


II. DIRECT KRAMERS-KRONIG ANALYSIS 
OF MAGNETO-REFLECTIVITY AND KERR 
ROTATION 


We start by considering the magneto-reflectivity and 
Kerr rotation spectra of a bulk sample at normal inci¬ 
dence. The complex magneto-reflectivity tensor has the 
same symmetry as the MO tensor 


r xx {u) r xy (u) 
-r xy {uj) r xx (u ;) 


( 6 ) 


Its eigenvalues are given accordingly by: 

r± (w) = r xx (uj) ± ir xy (uj) = \JR± (7) 
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In the local limit they are related to the momentum- The sought phases can now be obtained by applying the 

independent dielectric function via the Fresnel equations: KK transformation to R+(w) and 


r± M 


i - \A±M 

i + \A± M 


( 8 ) 


0±(oj) 



In \JR±{x) 
x — UJ 


dx + 7r, 


(15) 


The magneto-reflectivity R(u> ), unlike the zero-field re¬ 
flectivity, depends on the polarization state of the inci¬ 
dent light even if the sample is optically isotropic. For a 
circular polarization, it takes one of the values R±(uj) de¬ 
pending on the handedness. For the linear polarization, 
which can be represented as the sum of the two oppo¬ 
sitely handed circular waves with equal amplitudes, it is 
given by the average value: 


R+(u}) + R-(u}) 

R (u) = -j-' 


( 9 ) 


which can be reduced, using the parity relation 
R±(—u) = Rz p(w), to: 

0±{u)=9(ui)tQk(u), (16) 


where 


= ‘M-O + M-O 



In R(x) 
x 2 — 0 J 2 


dx+TT (17) 


and 


For a general (elliptical) polarization, the reflectivity 
is an unequally weighted superposition of i?_(w) and 
R + (uj). Notably, many spectrometers, especially of the 
Fourier-transform type, generate light with an ill-defined 
and frequency-dependent polarization. This imposes 
placing an extra polarizing element in a magneto-optical 
experiment in order to have a well determined polariza¬ 
tion. Below we shall consider the case where the incident 
light is linearly polarized and use Equation ([9 ). 

It is obvious that the standard KK method based on 
Equation ([5]) does not apply since the reflectivity phases 
0 + (uj) and 0_(w) are different. Moreover, it does not use 
the essential information coming from the Kerr angle: 




— 0 + (co) 


( 10 ) 


Here we prove mathematically that the complex re¬ 
flectivities r±{u>) can be model-independently extracted 
from R(uj) and 9k(u) if these function are known on the 
entire real-frequency axis. We first introduce an auxiliary 
complex function: 

/M = h 1 (j^j) = ^ lnc M + 2iO K {u}), (11) 


where 


c(w) 


R-M 

R+ M 


( 12 ) 


is the circular-dichroism ratio. As the imaginary part 
of /(w) is the Kerr angle known from experiment, 
the circular dichroism can be computed via the KK 

transformation:® 


c(w) = exp 



9k(x) 

x 2 — UJ 2 



(13) 


Combining this result with Equation © we can extract 
separately the reflectivity for each handedness: 


i?h = = rr^y fi(w) - (18) 


Equations (13)-(U8| constitute a self contained recipe 


to extract the complex MO reflectivity for each hand¬ 
edness and Equation ([8]) allows computing the complex 
dielectric functions e±(w). An important implication of 
this result is that spectra of i?(w) and 9k{u) measured 
in a broad spectral range and supplemented with proper 
extrapolations, bear sufficient information for complete 
MO spectroscopy. 

It is worth noting that the described algorithm relies 
on the non-trivial assumption that the functions In r± (w) 
are analytical, i.e. they do not exhibit poles or other 
singularities in the upper complex semiplane > 0. 
While this holds for normal reflectance from a bulk sam¬ 
ple, the analyticity can break down, for example, for a 
partially transparent sample. In this case, the variational 
KK constrained analysid® should instead be applied. In 
the next section we describe the application of the varia¬ 
tion approach to a much broader set of magneto-optical 
experiments. 


III. VARIATIONAL KRAMERS-KRONIG 
ANALYSIS OF MAGNETO-OPTICAL SPECTRA 

First we remind the basic idea of the KK constrained 
variational analysis in zero field®. The measured spectra 
such as reflectivity, transmission or any combination of 
them are initially fitted with a KK consistent and sum- 
rule satisfying model dielectric function. Most typically, 
this is a sum of a limited number of Drude-Lorentz oscil¬ 
lators: 


^mod(^) — ^oo + ^ ) 


Up,k 

w o,fc - w 2 - iyfcw ’ 


(19) 


R+ M 


n , ' 2c(w)R(u) 

1 + c(u>) ’ 1 + c(w) 


(14) 


where each Lorentzian is parametrized by the oscillator 
frequency wo,fc (equal to zero for a Drude component), 
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the linewidth 7 k and the plasma frequency uj Pt k- The 
parameter is the contribution from the high-energy 
optical transitions. Once a reasonable fit is achieved, 
these parameters are fixed, and a so-called variational 
dielectric function (VDF) is added to the model: 


e(w) — e mo d(w) T c var (uj). (20) 

The VDF is defined on a dense mesh of the experimental 
frequency points 07 (* = 1, ..,-/V), which usually coincide 
with the experimental datapoints or two times less dense 
(i.e. every second datapoint). It is constructed as a sum 
of elementary KK-consistent functions e^(ui): 


with the following structure: 

N-l 

£var±(w)= [ j4 i± £ fM+^F e f*(~ W ) 

i =2 


(26) 


that satisfies the physical relation e±(— uj) = e^(w). 
One can see that now there are two adjustable param¬ 
eters for every mesh frequency: Ai + = e val +( 07 ) and 
Ai_ = e va r-(wi). By refitting the data we obtain model- 
independently the complex valued e+(w) and e_(w) and 
thus achieve the goal of MOKKA. 


IV. EXPERIMENTAL PROCEDURE 


N -1 


i=2 


M = Ai e f M + e f*( _w ) 


( 21 ) 


where the coefficients A,; = e var (uJi) are regarded as free 
parameters. The imaginary part of ef'(w) has a triangu¬ 
lar shape: 


Sef (w) = 


UJi—LJi — i 
OJi+l — U) 


, 1 < UJ < Ui 

, LOi < UJ < Cj i+1 ( 22 ) 


Ui+1— UJi 

0 , otherwise 


while its real part is obtained by the KK transformation: 


1 f c 

M = -p 
n J-, 


^±dx =I 


g{u> - UJi- 1) 


UJi — UJi -1 


{ujj+i - 07-i)ff(u; - uij) g(o; - 07+1) 

(UJi - Wi_l)(w i+ i - Wj) Wj+l - UJi 


(23) 


where < 7 ( 0 :) = xln|x|. The second term to Equation (21) 
serves to satisfy the physical condition e(— uj) = 

The entire collection of data is finally refitted by adjust¬ 
ing all Ai values independently. As a rule, the fit quality 
is nearly perfect since the number of parameters is equal 
or comparable to the number of datapoints. Since the 
variational method does not imply the analyticity of any 
other response function than e(ui) it applies to virtually 
any type of optical experiments. 

It appears that this approach can also be applied to 
magneto-optics. The Drude-Lorentz model is modified 
in the magnetic field as follows: 


£mod,±(w) — £oo± T ^ ( 


UJ 


p,k 


k w o,fc 


■ UJ 1 T uj c ^kUJ 


, (24) 


where a cyclotron frequency uj c ^ is introduced for each 
Lorentz term. At first one has to fit the entire set of MO 


data as close as possible using Equation (24) or any other 


KK-consistent and sum-rule compliant function. Then 
the model parameters should be fixed and a magneto- 
VDF should be added: 


£±(^) — ^mod±(w) “f“ £ V ar± (^0 


(25) 


In this section we shall describe an experiment 
to measure broadband magneto-reflectivity (magneto¬ 
transmission) and the MO rotation spectra that can be 
treated by MOKKA. Although the measurement proto¬ 
cols described below were optimized for Fourier trans¬ 
form infrared (FTIR) spectroscopy, they can also be ap¬ 
plied, with some modifications, to other types of broad¬ 
band spectrometers. We assume that the experimental 
setup allows the inversion of the magnetic field direction. 
Hereafter until the end of the Section we shall refer to the 
reflectance and the Kerr rotation, although exactly the 
same experimental procedure applies to the transmission 
and the Faraday rotation. 

As it was emphasized in Section [TTJ in magneto-optics 
the polarization state of the incident light should be well 
defined even if the sample is isotropic in zero field. Since 
the spectrometer output may be uncontrollably polar¬ 
ized, a polarizer should be put before the sample. Ad¬ 
ditionally, the optical path between the sample and the 
detector may also introduce some depolarization effects, 
in which case placing a second polarizer (often called an¬ 
alyzer) after the sample is also required. Therefore we 
shall consider a setup, where both polarizers are present, 
and which allows measuring the reflectivity and the Kerr 
rotation within the same experiment. This configura¬ 
tion is similar to standard two-polarizer ellipsometrj®, 
except that here we consider normal incidence, while el- 
lipsometry is typically done at large illumination angles. 
We shall assume that the polarizer is fixed and the ana¬ 
lyzer is rotated, although the same result is valid for the 
case where the analyzer is fixed and the polarizer turns. 
Since the sample is isotropic, we can set the polarizer an¬ 
gle to zero without loss of generality. As it is shown in 
Appendix A, the intensity measured by the detector as a 
function of the analyzer angle A is given by: 

I A {uj, B) = C a (uj, B) 

R{uj , B) + R(uj , B ) cos[2A - 2 0 k (uj, B)] /0 ^ 

x 2 > lA'I 

where the functions R(co,B), R(uj,B) and 6k{uj,B ) are 
defined in Section [TT| and the prefactor Ca(uj, B) is deter¬ 
mined by all elements in the optical path from the source 
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to the detector apart from the sample. It may depend on 
the analyzer angle due to depolarization effects in the de¬ 
tection channel and on the field because of the changes in 
detector sensitivity and all mechanical movements that 
it may induce. 

If we replace the sample by a reference with R = R = 1 
and 9k = 0 then the detector intensity will be given by 
the Malus law: 


I A ,ref{u, B ) = Ca(u, B) cos 2 A. (28) 


Thus the reference measurement allows us to determine 
C A (ui, 33)P2. Since the choice of the best reference proce¬ 
dure is setup- and sample-dependent, we do not discuss 
this issue her^H. 


Below we discuss two measurement protocols, both 
based on Equation (27) but adapted for different experi¬ 
mental situations. 


A. Complete protocol 


As it is common in polarimetric experiments, the most 
rigorous approach involves measuring the raw spectra 
I a {w,B) for a large number of analyzer angles spread 
equidistantly between 0° and 360°. The same should 
be done for the reference, in order to obtain the prefac¬ 
tor C a (uj, B ) as discussed above. For each frequency and 
each field value one can fit the function F(A) = I A /C A by 
the formula a + (3 cos(2A — 2 A m ) with the fitting param¬ 
eters a , /3 and A rn . The reflectivity can be determined 
by the relation R(uj,B) = 2 a(uj,B). The parameter A m 
formally corresponds to the Kerr angle. However, using 
the relation 6bf(w, 33) = A m (co,B) may result in a con¬ 
siderable systematic error due to polarizer imperfections. 
This problem is easily removed by subtracting the zero- 
field value: 9k{<jj,B) = A m (uj,B) — A m (uj, 0). 

Our experience shows that measuring at the two po¬ 
larities of the magnetic field improves both the system¬ 
atic accuracy and the signal-to-noise ratio. Based on the 
parity relations 3?(w, —33) = i?(w,33) and 9k(u,—B) = 
—9k(w,B) the following equations can be derived: 

R(u, B) = a(ui, B) + a(u, — B), (29) 

9 k {uj, B) = A ™(u,B)~A m (u,~B) . (3Q) 

that we actually recommend to use. 

The value of /3/a = R(u>)/R(u>) = 2y/c(u>)/(l + c(w)) 
can in principle be utilized to determine the circular 
dichroism ratio c(w). However, this only works well if 
the dicroism, and therefore the ellipticity of the reflected 
light, is strong enough, which is often not the casd^. 
Fortunately, the implementation of MOKKA does not 
require c(w) since it is computed via a KK transforma¬ 
tion of 0k(w). 


B. Fast protocol 


The measurement scheme described above is rigorous 
but quite slow as it involves measuring the raw intensity 
spectra at dozens of analyzer angles for the sample and 
the reference. Fortunately, in many situations one can 
use an alternative protocol, which requires much fewer 
spectra to be measured. This method is accurate if the 
rotation angle is small (for example 9k < 0.2 rad), so 
that sin 9k ~ 9k and cos 9k ~ 1. Since for small rota¬ 
tions the dichroism is also weak, \c — 1 | -C 1 , the ratio 
R{w)/R(w) = 1 + 0((c — l) 2 ) and hence we can ignore 
the small difference between R(ui) and R(lo). 

The fast method hinges on the fact that the intensity 
is most sensitive to the MO rotation when the polarizers 
are oriented at +45° or —45° with respect to each other 
and is almost insensitive to it for A = 0, i.e. when they 
are parallel. Within the approximations made above, 
Equation (271 gives the following detector intensity for 
these analyzer positions: 

3±45( w ,33) ~ C±4$(uj,B)R(uj,B) 
Io(u,B)nC 0 (u,B)R(uj,B). 
and for the reference: 

C±45(w, B) 


1±2 9 k (u,B) 


I± 45 ,re f{u,B) 

(kb 33) 


Co(w, B). 


,(31) 

(32) 

(33) 

(34) 


An obvious way to determine the reflectivity is to use 
the measurement in the parallel polarizer configuration: 


R{oj,B) 


Iq(u>,B) 

hi.ref (k>, B) 


(35) 


which has the additional advantage that the optical signal 
is at maximum. 

In order to obtain the magneto-optical rotation, we 
first we assume that Ca(w, +B) = Ca(w, — B ). Then we 
have: 


P±45(u, B) = ± 4 °j; ’ + \ ~ 1 ±46 k (u,B), (36) 
3±45(w,~B) 

where again we used that 9k(w,—B) = — 9k(oj,B ). This 
suggests that the rotation can be obtained based on the 
measurement for only one of the angles +45° or —45°, 
without reference: 


9k (w, B) ss ± 


P±45(v,B) - 1 

4 


(37) 


In this way the rotation can be measured even in exper¬ 
imental setups, where in situ rotation of the analyzer is 
not possible as one can simply place the sample between 
two fixed polarizers mounted at 45° with respect to each 
other. 

If Ca(u>, +B) is not equal to C a (uj, —33) then the pro¬ 
cedure should be modified. The most obvious modifica¬ 
tion is to take the same measurement with a reference: 


P±45,re MB) 


3±45 ,re f(v,+B) 
3±45,re/(kb 33) 


(38) 
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and to use the corrected formula: 


e K (u,B)^± 


P±45(w,-B) 


1 


4 P±45,re /+,£) 


- 1 


(39) 


If the analyzer can be rotated then a more accurate 
result can be obtained by combining measurements at 
+45° and —45°. One can use the formula: 


Ok(u,B) 


P+45+, B) - p_ 45 (u,B) 


(40) 


if Ca(w,+B) = Ca{w,—B) or 


9 k (oj,B)^1 


P+45+, -B) 


p 45+,B) 


P+45,re f((jJ,B) P—A5,ref + , B) J 


(41) 


otherwise. Equations (WOfc and (411 are less sensitive 


to polarizer imperfections than Equations (371 and (39). 


Additionally, one can show that Equation (40) provides 


a more accurate result than Equation (37) when a weak 


field asymmetry of Ca+, B) is present. This allows a 
reference-free measurement of the MO rotation without 
imposing stringent conditions on the stability of the op¬ 
tical elements with respect to the field. 


V. APPLICATION EXAMPLES 
A. Bismuth 






FIG. 1. Far-infrared optical and magneto-optical infrared 
spectra of bismuth at 10 K for the polarization perpendicular 
to the trigonal ( z ) axis. Panels (a) and (b) show the zero- 
field reflectivity and optical conductivity respectively. Panel 
(c) presents magneto-reflectivity (in blue) and Kerr rotation 
(in red) at 7 T. (d) The real (absorptive) part of the optical 
conductivity for the two circular polarizations obtained using 
MOKKA. The inset of panel (b) depicts the electron and hole 
Fermi pockets (adapted From Ref.' 18 ). 


Bismuth is a canonical semimetal with a low car¬ 
rier density (~10 17 cm -3 ), a long mean-free path and a 
small cyclotron mass. These properties, combined with 
a strong spin-orbit interaction, give rise to spectacular 
magneto-optical effectipZHSil While some classical phe¬ 
nomena, such as quantum oscillation d 32 * 33 ! were first dis¬ 
covered in bismuth, a number of new phenomena were 
recently found in it and in related compounds, such as a 
transition to a 3D topological insulating state upon Sb 
substitutiorPS, a valley-ferromagnetic statJ^, a valley- 
nematic Fermi liquid statd^ and avoided Lifshitz-type 
semimetal-semiconductor transition^. Therefore this 
material makes an excellent case for applying MOKKA. 

Bismuth has a rhombohedral (or trigonal) crystal lat¬ 
tice and therefore it is optically anisotropic. However, it 
is optically isotropic within the plane, which is perpen¬ 
dicular to the trigonal axis and therefore the described 
technique applies when the propagation of light and the 
magnetic held are along the trigonal axis. Fortunately, 
the material cleaves easily (in liquid nitrogen) in the same 
plane, making it possible to prepare a high-quality atom¬ 
ically flat surface of the needed orientation. 

Infrared magneto-reflectivity and Kerr angle spectra at 
nearly normal incidence (~ 8°) were measured between 
20 and 650 cm -1 with the resolution of 1 cm -1 at 10 
K in a split-coil superconducting magnet attached to an 
FTIR spectrometer. Both protocols described in Section 
|IV| were used for comparison and they gave almost iden¬ 
tical results. A Hg light source and a liquid He cooled 


bolometric detector were used. A crystal of 99.9999% 
pure Bi had a surface of 3 x 5 mm 2 and a thickness of 1 
mm. Free standing wire-grid gold polarizers were used. 
The magnetic field was applied along the trigonal axis 
and therefore almost parallel to the propagation of light. 
The noise in the reflectivity spectra was less that 0.5 %. 
The accuracy of the Kerr angle was better than 1 mrad. 

The zero-field reflectivity, shown in Figure [1^, features 
a sharp plasma edge at around 150 cm -1 due to highly 
mobile charge carrie rs and it is in good agreement with 
previous workP^Ml The optical conductivity obtained 
by the usual KK analysis, is presented in Figure [T]). It 
has a low-energy Drude peak and a linearly growing con¬ 
tribution from low-lying interband transitions. 

The magneto-reflectivity and the Kerr angle spectra at 
7 T are shown in FigureJTj-. One can note that the plasma 
edge is significantly blue-shifted by the field, which is 
a magnetoplasma effeclP- 1 . The reflectivity also shows a 
strong decrease below 100 cm -1 , which can be attributed 
to the enhancement of the penetration depth and there¬ 
fore an enhanced optical absorption below the cyclotron 
frequency. The Kerr rotation has a strong spectral struc¬ 
ture clearly associated with the plasma edge at 200 cm -1 
and a weaker feature at 100 cm -1 , close to the frequency 
below which the reflectivity starts to decrease. Remark¬ 
ably, the Kerr angle is essentially zero above 250 cm -1 . 

Figure [ljl presents the optical conductivity a± +) for 
the left- and right-handed circular polarized light ob¬ 
tained by the variational magneto-optical KK analysis, 
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where a Drude-Lorentz model with only a few terms was 
used as a starting point. Applying the direct integration 
procedure described in Section [IT] provided essentially the 
same conductivity curves. The most prominent spectral 
structure in both polarizations is a absorption band at 
80-110 cm -1 , which consists of several peaks. The spec¬ 
tral shape of this absorption is quite different for the two 
light helicities. The Fermi surface in bismuth is highly 
anisotropic. It features one hole pocket and three elec¬ 
tron pockets with a Dirac-like band dispersion as shown 
in the inset of Figure [I] >. The presence of multiple cy¬ 
clotron energies is obviously related to this complex band 
structure. 

The physical interpretation of this data is beyond the 
scope of this paper and will be presented elsewhere®! 
Here we only note that MOKKA allows us to obtain the 
true shape of the MO conductivity and to determine the 
absolute spectral weight of the cyclotron resonance peaks 
for each handedness separately in a model-independent 
fashion. 


B. Graphite 

Graphite is another textbook semimetal, where several 
physical phenomena were discovered, and which is addi¬ 
tionally of a tremendous practical importance. In the last 
decade it again attracted much interest in the research 
community because of the exfoliation of graphene®!. 
Graphite is composed of weakly coupled honeycomb car¬ 
bon layers with predominantly Bernal stacking and is 
therefore extremely anisotropic. The tiny Fermi surface 
consists of trigonally warped electron and hole pockets 
stretched along the K-H corner of the hexagonal Brillouin 
zone, as shown in the inset of Figure [2 Jd. The very pe¬ 
culiar optical and magneto-optical properties of graphite 
were widely studied 44 ®!. Recently, we measured infrared 
magneto-reflectivity and Kerr rotation spectra and fitted 
them using a tight-binding model HamiltoniarP^. Here 
we shall demonstrate a model-independent extraction of 
handedness-resolved magneto-optical conductivity from 
the same data. 

Figures [5ji and[2]D present the far-infrared region of the 
zero-field reflectivity and the KK-derived optical conduc¬ 
tivity of highly-ordered pyrolytic graphite (HOPG) at 10 
kP^I. The reflectivity is growing monotonously towards 
unity at zero frequency and the optical conductivity is 
essentially constant, apart from the lowest frequencies, 
where a small Drude peak is observed. Above 200 cm -1 
the conductivity is close to do = e 2 /4 h per layer (« 1820 
H _1 cm _1 ), which is the theoretically expected conduc¬ 
tivity of monolayer graphene due to optical transitions 
between conical valence and conduction bands®!. 

Figure [2]:’ shows the experimental results at 7 T. The 
reflectivity modifies dramatically in field showing intense 
spectral structures owing to the formation of Landau lev¬ 
els. The Kerr rotation also shows a strong and rather 
nontrivial spectral dependence. We successfully modeled 



FIG. 2. Far-infrared optical and magneto-optical spectra 
of graphite at 10 K. Panels (a) and (b) show the zero-field 
reflectivity and optical conductivity respectively. Panel (c) 
presents magneto-reflectivity (in blue) and Kerr rotation (in 
red) at 7 T. (d) The real (absorptive) part of the optical 
conductivity for the two circular polarizations obtained using 
MOKKA. The dashed curves in Panels (c) and (d) correspond 
to the tight-binding Hamiltonian fitfP. The inset depicts the 
Fermi surface (electrons pockets in blue, hole pockets in red), 
adapted From Ref.®). 


these data in fields from 0 to 7 T using a Slonczewski- 
Weiss-McClure Hamiltonian with Landau-level quanti¬ 
zation and the Kubo formula as described in Ref. [49] 
The fits at 7 T are shown as dashed lines. Although 
very good, they show some deviations from the experi¬ 
ment and therefore do not provide the precise shape of 
the optical conductivity. Here we use these physically 
meaningful fits, instead of the phenomenological Drude- 
Lorentz model, with fixed parameters as a starting point 
and apply a magneto-VDF as discussed in Section |III| 
This allows us to achieve a nearly perfect fit (not shown) 
and calculate the optical conductivity for the two circu¬ 
lar polarizations shown in Figure [2ji. One can see that 
the true optical conductivity shows much sharper peaks 
than the model. This is likely due to the fact that the 
Kubo formula used in Ref. US] was assuming a frequency- 
independent scattering rate, thus strongly overestimating 
the optical width of low-energy LL transitions. 


In contrast to the case of bismuth, the two conduc¬ 
tivity curves of graphite are markedly different. This is 
explained by non-trivial magneto-optical selection rules 
in graphite that are dominated by the trigonal warping, 
as discussed in details in Refill 















c. 


Multilayer graphene on SiC 


As a third example, we show an extraction of 
handedness-resolved magneto-optical conductivity of an 
ultrathin film on a MO-inactive substrate by applying 
MOKKA to a combination of magneto-transmission and 
Faraday rotation. Figures [3^ and [3}} present far-infrared 
spectra of T(B)/T( 0) and dp of multilayer graphene 
grown on the C-face of 6H-SiC 51 for several fields be¬ 
tween 0.25 and 4 T at 5 K. The Fabry-Perot oscillations 
in the substrate were intentionally suppressed by reduc¬ 
ing the spectral resolution to 5 cm^ 1 . 

The two-dimensional (sheet) optical conductivities 
ct±(uj) of graphene (Figures |3j: and[3ji)) were calculated 
model-independently using the variational MOKKA in¬ 
troduced in section m and the Fresnel formulas for the 
case of incoherent internal substrate reflections as de¬ 
tailed in Appendix B. 

The transmission features several absorption dips 
(marked by arrows for B= 0.25 T) corresponding to opti¬ 
cal transitions between Land au levels, in agreement with 
previous measurements®!®!. The transition energies are 
proportional to \[B and match almost perfectly the ex¬ 
pected values in monolayer graphene, where the LL en¬ 
ergies are given by the formula: 

E n = sign(?r) 2heVp\nB \. (42) 




D + 

<D 

cn 



Here vp ~ 10 6 m/s is the Fermi velocity and n = 
0, ±1, ±2,.. is the LL quantum number. It was argued®! 
that the monolayer-like dependence is due to a twisted 
stac king of the layers which decouples them from each 
other®! much more efficiently than in the Bernal-stacked 
graphite. The absolute spectral weight of the LL absorp¬ 
tion is in fact much smaller than the expected value ac¬ 
cording to the number of layer (about 20 in this sample) 
for reasons that are not fully understood at the moment. 
In addition to quantum Landau level transitions, there is 
a structure at low energies (shown by the asterisk for 4 
T) stemming from the cyclotron resonance of free carriers 
in strongly doped layers close to the substrate. 

While the LL transitions —n —> n + 1 and —n — 1 — > n 
have the same energy, they are excited by circular polar¬ 
izations of light®! as sketched in the insets of Figures [3jc 
and|3j:l). In a simple magneto-transmission experiment 
they unavoidably overlap. The advantage of MOKKA 
over is that it allows an unambiguous separation of transi¬ 
tion with different handedness. In the presented case, one 
can clearly see that the shape of the peaks with the cor¬ 
responding energy is quite different in cr + (u;) and er_(u;). 
For example, the peaks in o+(w) are sharper. One can 
also notice that at low fields (0.25 and 0.5 T) the inten¬ 
sity of the 0 —> 1 peak is significantly larger than that 
of the — 1 —> 0 peak. A discussion of this electron-hole 
asymmetry was given in RefP. 


FIG. 3. Far infrared magneto-transmission (a) and Faraday 
rotation (b) of multilayer graphene on the Si-face of 6H-SiC 
at 5 K at several values of perpendicular magnetic field be¬ 
tween 0.25 T and 4 T. Arrows indicate the LL transitions 
from neutral layers, the asterisk shows the Drude contribution 
from highly doped layers. Panels (c) and (d) show the two- 
dimensional optical conductivity normalized by <to = e 2 /4 ft 
for circular light of each handedness obtained using varia¬ 
tional MOKKA. The corresponding Landau level transitions 
are sketched in the insets. In all graphs, the curves are shifted 
vertically for clarity, except at 0.25 T. The zero levels for the 
Faraday rotation curves are shown by the dashed lines. 


VI. SUMMARY 

The goal of this paper is to introduce a technically sim¬ 
ple solution for a complete-tensor broadband magneto¬ 
optical spectroscopy. The pathway that we propose is to 
perform an ellipsometry-like measurement at normal in¬ 
cidence with two polarizers and to supplement it with a 
properly generalized Kramers-Kronig analysis. We have 
proven mathematically that in the case of a bulk sam¬ 
ple a combination of broadband magneto-reflectivity and 
Kerr angle that can be measured using this simple tech¬ 
nique, allow a unique and model-independent extraction 
of handedness resolved dielectric function. Furthermore, 
we extended the variational Kramers-Kronig analysis 21 
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to magneto-optical spectra, which can be applied in a 
variety of experimental situations, including transmis¬ 
sion and Faraday rotation measurements on thin films 
deposited on a substrate as well as reflectivity on semi¬ 
transparent samples. This method can be also regarded 
as an optical-Hall type technique since it provides access 
both to the longitudinal and the Hall components of the 
optical conductivity. 

As the examples of bismuth, graphite and epitaxial 
graphene on SiC have clearly shown, handedness resolved 
MO spectroscopy is especially useful for the systems, 
where both electrons and holes are present. In contrast 
to the often used approach, where the measured spectra 
(including Kerr or Faraday angles) are only compared 
to model curveiP^ESdll^ our technique allows extracting 
the true shape of the MO dielectric function and opti¬ 
cal conductivity. This is important, where such a shape 
is difficult to reproduce theoretically, for example, due 
to a frequency-dependent scattering rate. Although all 
spectra presented in this paper were measured in the far- 
infrared region, the same technique can be obviously ap¬ 
plied to any other optical range. 

Further extensions of this method to more complicated 
situations can be envisaged. First, in the case, where 
the magnetic and magneto-electric susceptibility play a 
role one would need to increase the number of indepen¬ 
dent measurements. For example, measuring magneto¬ 
transmission and Faraday rotation for several films of 
different thickness and properly modifying the Fresnel 
equations may allow an independent extraction of e± (uj) 
and h±(uj). Second, if non-local effects are essential 
then one should introduce the momentum dependence 
to the dielectric function and replace the Fresnel equa¬ 
tions with more general electromagnetic relations. Re¬ 
markably, even in this case the complex reflectivity for 
each handedness can be extracted using the formalism 
presented in Section [Hj which does not rely on the Fres¬ 
nel equations. Third, if a sample is anisotropic in zero 
field and the MO tensors e(w) and r(ui) no longer sat¬ 
isfy the symmetry of Equations (|2| and © the measure¬ 
ment technique has to be modified to involve rotating of 
both polarizers. As a final remark, the time-domain po- 
larimetric spectroscopj^MH offers further opportunities 
when combined with the method described in this pa¬ 
per, as it allows measuring more independent quantities 
simultaneously. 


Acknowledgments 

We thank M. Tran, N. Ubrig, D. van der Marel, J.-M. 
Poumirol and A.Akrap (University of Geneva) for valu¬ 
able discussions. We are grateful to C. Uher (University 
of Michigan) and T. Seyller (University of Chemnitz) for 
providing us high-quality samples of bismuth and epi¬ 
taxial graphene on SiC. This work was supported by 
the Swiss National Science Foundation and EU Graphene 
Flagship (contract no. CNECT-ICT-604391). 


APPENDIX A 


Here we use the Jones-matrix formalism to derive 
Equation (271. For simplicity, we shall refer to the re¬ 
flection experiment, although exactly the same deriva¬ 
tion applies to transmission. Below we consider the case 
where the first polarizer is fixed and the second polarizer 
(analyzer) is rotating. However, the same result applies 
when the first polarizer is rotated and the second one 
is fixed. The transformation of the polarization of light 
along the optical path is sketched in Figure [3] 



FIG. 4. Transformations of the polarization of light in the 
two-polarizer optical setup used for the magneto-optical ex¬ 
periment proposed in this paper. 

Without loss of generality, we assume that the polar¬ 
izer is aligned along the x axis, so that its Jones matrix 
is 


P = (J S) <*» 

where tp is the polarization transmission for the wanted 
polarization. This results in the following Jones vector 
of the radiation incident to the sample: 


Ein — PE source — tpE °’ 

where Eq is the x component of the amplitude of light 
leaving the spectrometer. Using the Jones matrix of the 
sample 



we obtain the Jones vector of the reflected light: 


Eref — SE^ 



(46) 


which in general corresponds to an elliptical polarization. 
The Jones matrix of the analyzer set at the angle A with 
respect to the x axis is: 


/ cos 2 A 
^cos A sin A 


cos A sin A\ 
sin 2 A ) ’ 


( 47 ) 
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where t A is the analyzer transmission. After the analyzer 
the following polarization state is achieved: 

Rout = AE re f 

sin a) ( rxx cos ^ “ rx v sin ^) tAtpE 0 , (48) 
The signal at the detector is finally given by 



I a (lo) = C E out = C A \r xx cos A-r xy suiA\ (49) 


where £ is the cumulated transmission coefficient of all 
optical elements between the analyzer and the detector 
and C A = (\t A tpE 0 \ 2 . Note that ( and C A may depend 
on A if some of these elements are polarization sensitive. 
In the following, we take into account that 


y /R + e~ ieK + y/Rle i8K 


^ XV — 


\/R+e 


— i6 b 


- s/RZe i9K 


2 i 


,i(9++6 K ) 

(50) 

,i(0 + +0Rr) 

(51) 


as it follows from Equation (J7|. Substituting this into 
Equation (491 yields the dependence of the detected sig¬ 
nal on the analyzer angle: 

(y/R+e~ i0K + ^Rle WK ) cos A 


T - Ca 


+i( v A^ e - ie K _ y /RZe if>K ) sin A (52) 


After simple algebra this reduces to: 

R + R cos(2A — 26 k) fe . Q ^ 

I a = C A ---■ (53) 

This result has a straightforward geometrical interpre¬ 
tation. The detector signal is minimized when the ana¬ 
lyzer is orthogonal to the main axis of the ellipse drawn 
by the electric-field vector of light reflected from the sam¬ 
ple (Figure [4]). On the other hand, the orientation of this 
axis is given precisely by the Kerr angle. 


T± 




Zpcr± 
n s + 1 


-2 


1 f n s + 1 + Z 0 cr + \ 

2 arg U + l + ^J 


(54) 

(55) 


where Zq is the vacuum impedance. In this situation the 
effect of the substrate as compared to a free standing 
film is to reduce the effective value of the film conductiv¬ 
ity by (n s + l)/2 due to a screening of the electric field 
experienced by the film. 

If multiple reflections contribute significantly to the 
signal then the effect of the substrate is more complex 
and is dependent on the degree of their mutual phase 
coherence. Decoherence can be caused by wedging of 
the substrate, a distribution of the angles of incidence, 
a finite resolution of the spectrometer and other factors. 
While it is in general quite hard to reproduce this effect 
quantitatively, the analytical result can be easily found 
in the two limiting cases. If the internal reflections are 
fully coherent then 


Tl = TV x 


1 — p 2 


r\c _ r \0 

Up — Up 


i 


i -pq± 

1 - pq+ \ 


arg 


1 -pq-J ' 


where 


P = 

q± = 


n s - 1 

—rr exp ' 

n s + 1 
n s — 1 — Z 0 cr± 
n s + 1 + Zq(t± 


{i“n s d s }, 


exp 


{i- c n s d s ). 


(56) 

(57) 


(58) 

(59) 


and the Fabry-Perot effect is present both in the trans¬ 
mission and the Faraday rotation—. In the opposite case 
of the full incoherence®* one obtains: 


APPENDIX B 

For the reference purpose, here we present equations 
needed to calculate the magneto-transmission T = (T + + 
T_)/2 and the Faraday rotation Op of a MO-active film 
with the two-dimensional optical conductivity cr± de¬ 
posited on a MO-inactive substrate with the complex 
refractive index n s and thickness d s . The result de¬ 
pends strongly on whether or not the multiple internal 
reflections inside the substrate are taken into account. 
If they are negligible, which is, for example, the case 
for a strongly absorbing substrate, then the transmission 
normalized to the bare identical substrate and Faraday 
rotation are given by: 


Tl= Tlx 


1 - \p\ 4 

1 - bid 


12 ’ 


9 i F = 0° F + -axg{l- \p\ 2 q*_q+). 


(60) 

(61) 


Interestingly, at selected frequencies the Fabry-Perot 
effect can increase the transmission and the Faraday rota¬ 
tion simultaneously as compared to the incoherent case, 
which can be important for MO application^®!. On the 
other hand, applying MOKKA is easier if the Fabry- 
Perot oscillations are fully suppressed. 
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